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We calculate the correlation functions of 0^^ qq, ss and glueball in the QCD sum rules and 
obtain the mass matrix where non-diagonal terms are determined by the cross correlations 
among the three states. Diagonalizing the mass matrix and identifying the eigenstates as 
the physical O^"*" scalar mesons, we can determine the mixing. Concretely, our calculations 



< 

QQ , determine the fractions of qq, ss and glueball in the physical states /o(1370), /o(1500) and 



/o(1710), the results are consistent with that gained by the phenomenological research. 

Ph. 

Dh' I. INTRODUCTION 

D 

Existence of glueball is a long-standing puzzle in the QCD theory. Searching for it becomes the 

(N : 

■ most challengeable task for high energy physics society. The QCD theory predicts its existence and 
OO ■ r 

the lattice QCD almost determines the mass spectra of glueballs with various quantum numbers[ll - 
9I. It is believed that the mass of the lighter glueballs should be at around 1 to 2 GeV. But where 



(N 

O ■ are they, can we pin down them? Several bound states near 2 GeV have been found in recent 
experiments |lo| . People believe that the number of these states indeed exceeds that predicted by 
- ^ the simple symmetry analysis. One natural explanation is that there exist exotic states and the 

^ '. newly observed resonances are either such exotic states, glueball, hybrid and multi-quark states, 
or their mixtures. In fact, none of the resonances which are newly observed at BES and BELLE 
can be identified as glueballs, so that one is tempted to conclude that glueballs mix with the regular 
quark states. The lattice and other model-dependent calculations all predict the mass of the 0+"*" 



glueball falling within the range of about 1.7GeV[|lll- ll3|l . Meanwhile the mass of the state made 



of pure light quarks qq, where, the q refers to u, d and s quark, is also near L3 ~ 1.7GeV[|l4l-ll6h. 



therefore it is very possible that the scalar glueball and the quark states mix to constitute physical 
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states. The observed resonances /o(1370), /o(1500) and /o(1710) which have masses close to 
[1.3 ~ 1.7]GeV, can really be such mixtures. 

For this mixing, it implies that the scalar glueball does not independently exist as a physi- 
cal state which people explore in experiment, but the three physical states: /o(1370), /o(1500) 
and /o(1710) possess glueball components. In fact, many authors have discussed the mixing of 
these three physical states [Il7l - l2 ill . Generally, this issue was discussed based on phenomenology, 
namely by fitting data of various reactions, the mixing parameters are fixed. It would be inter- 
esting to investigate this problem from a more fundamental theory. However, the energy scale 
for the mixing is low and the non-perturbative QCD effects may dominate, therefore the regular 
perturbative theory does not apply. By contrast, the QCD sum rules may be the bridge between 
perturbative quantum field theory and the non-perturbative phenomena i22|] . thus should be a rea- 



sonable approach for this research. Two groups have done the significant work [|23 



2511 . Narison 



et al's work fixed the mixing of the three states: /o(1370), /o(1500) and /o(1710) through the de- 
cays of the light-quark meson and the glueball. In their work, the masses of the scalar light-quark 
states and glueball are determined in the QCD sum rules and by using them to estimate the de- 
cay rates of the corresponding processes they fix the mixing parameters. By contrast, we assume 
that the scalar light-quark states | A^, S) and the glueball \G) are un-physical, therefore the masses 
independently determined in the QCD sum rules cannot be used to estimate the decay rates. In 
another work, Steele et al. predicted that the mixing states should involve mixing of /o(980) with 
the /o(1500) and /o(1710) in terms of the Gaussian QCD sum rule. Instead, in our work, we are 
going to investigate the mixing of the three states all near 2GeV in the QCD sum rules. 

The first step of our work is to define the currents for the un-physical states: glueball \G), light- 
quark states I A^) and 15) (N is for u, d quarks, and S is for s quark), then find their relations to the 
three physical states: |/i), I/2) and I/3) via a mixing matrix V. 

The work is organized as follows. After this introduction, we calculate the correlation functions 
in terms of the QCD sum rules, in Section III, we formulate the mixing matrix and show the 
relations between the unphyiscal states and the physical scalar mesons. In Section IV, we present 
our numerical results and the last section is devoted to our conclusion and discussion. 



II. THE CORRELATION FUNCTION 

In the scenario of the QCD sum rules, the correlation function n(g^) is defined as: 

nfe^,^./..e-(0|T{.W,.(0)}|0). (1) 
By the dispersion relation, at the hadron hand, the correlation function can be written as: 

n(,^, = i / (2) 

TT J s — q-^ 

After the Borel transformation and considering the quark-hadron duality, we obtain the "Moment" 
TZ as: 

1 

'Rk = - dssHmR{s)e-'^ (3) 

where r is the Boral parameter and sq is the threshold for the continuity. 
So in our work, the relevant correlation functions are defined as: 



where Jg(x) is 



n^^{q^)=i j rfxe^''^(0|T{Jq(x), Jq(0)}|0) 
n-(g2) = z j dxe^5-(0|T{Js(x),Js(0)}|0) 

W^{q^) = % J dxe^''"(0|r{Jg(x), Jg(0)}|0) (4) 
n^g(g2) = , J dxe^^^{0\T{Ux),J,{Q)}\Q) 
= , J da;e^5-(0|T{Js(x),Jg(0)}|0) 

Jg(x) = (5a) 



and J(i,s{x) is: 



Jq,s{x) = mq^sijq^s{x)i)q,s{x). (5b) 
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FIG. 1: The Feynman Diagrams for IlqgSjg :(a) perturbative part; (b-c) with quark condensates ; (d) with 
gluon condensates; (e) with quark-gluon condensates 



The "Moments" TZ are defined as: 
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(6) 



where 7^f can be found in Refs-OHH and 7^^''''' is given in Refs.|ll4|,|27l]. For the mixing 
current, we calculate the correlation functions and the "Moments" TZ^^'^^ are obtained from the 
Feynman Diagrams in Fig-[T1 

With the Operator Product Expansion (OPE), the correlation function Il'^'^'^^{q'^) is decomposed 



as: 



nqq.ss(g2) ^ Codo + Cs{qq,ss) + C,{a,G^)+C,{gsO,) 



(7) 



where Ci{i = 1, 3, 4, 5, ■ ■ ■ ) are the Wilson coefficients, and the operator Oq is the unit operator. 
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FIG. 2: The Counter temi for Co, where ■ is for the vertex correction 



In the fixed-point gauge [12 811 . we calculate the two-loop diagram in Fig-[TJa) and then we have: 



Co 



13a2 2 Q 
-m„ „ lor 



e 27r 



-77-0 ? 



log^ 



31og47r + 37b - ^ 



(8) 



where, = —q^. In Eq.®, we drop out the terms which are not proportional to Xo^Q'^ lv^\ 
because they do not contribute to the moment 1^^^ and disappear after the Borel transformation. 
The Feynman diagrams related to the counter terms are presented in Fig-l2l 



In the MS scheme, we find: 



(a) 



(9a) 



and 



(6) 



13a: 
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Eventually we have the coefficient Co at the two-loop order as: 



+ 



Co + Ct^ + Ci^^ 



Q^s 2 ^2 



3 13, 



(9b) 



(10) 



4 z^^ 2 ~ z^" 

which corresponds to the perturbative contribution to the moments. The other Wilson coefficients 
are calculated from Fig-fllb-e) as: 



Q,2 g2 
C3 = -47r— mq,s log -IT] 



a 



q,s 



a. 



-log:V + 3- 



z/^ 



TT 



(11) 



C5 = -^mnsOs. 



2 ""qjS^s- 
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With the correlation function Eq.©, the moment is: 

1 



where. 



— (1 -pi(sor))a[}'' - — ^ 



-IE + ^i(sor) + logsor + e 'V - 1 



-(1 - pi(sor))log 



■So 



hi 



;i - po{soT))mg{qq) 



+ 



Co' - Ci''(7i? + logrz/^ + ^i(sor)) 



(12a) 



q,s 



^^Qg 2 q,s _ 2 

4^3 ^g,s «i - "^9,' 
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a 



q,s 
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-m, 



vr 



q,s 



(12b) 



TT 



iq,s 



and pi,2 - (a;) and -E'i(x) are already given in [HQ 



It is noted that our result is different from that given in [|25|] . This is understood since different 
subtraction schemes are employed in the two works. 



in. EQUATIONS FOR MIXING MATRIX V 



We define the physical states as I/2) and I/2), whereas the un-physical states as | A^) = \qq), 
\S) = \ss) and \G). The mixing matrix connecting them is: 



1/2) 

[\h)J 



/t 



\s) 



(13) 



Vu V12 Vis 
V21 V22 V23 

According to the first approximation, it is assumed that |A^), \S) and \G) constitute a complete 
basis 11 1811 . but as a matter of fact, when the other resonances /o(1790) and /o(1812) were observed 



by the BES collaboration BO 



the aforementioned states [121 



3111. we suggested that the hybrids might join the game and mix with 



3211 ■ But it seems that one can first ignore the hybrids which might 



be heavier than the other three, and assume that the three physical mesons are only composed of 
the regular quark and glueball components. We will discuss this issue in the last section. So, the 
mixing matrix V transforms the flavor representation into the physical representation, i.e. the mass 
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representation, so it must be unitary, thus we have: 



11 ^ ^21 ^ ^-il 

V,\ + VI + VI = 1: (14) 



12 ''22 ''32 ~ 

13 + ^3 + ^33 



and the conditions are enforced 



^11^12 + ^2lV^22 + V^3lV32 = 0; 

V^iiV^i3 + ^2iV^23 + V^3iV^33 = 0; (15) 

^11^^13 + V2lV2-i + V^31^33 = 0. 

Next, we will build the equations to solve this mixing matrix V in terms of the QCD sum rules. 

In QCD sum rules, the integrand of the dispersion integral includes the imaginary part of the 
correlation function 11 at > and then one inserts a complete set of physical states of O"*""*" 



hadrons between the currents [|33ll . In the quark-hadron duality the lowest states' contributions 
dominate and the contributions of the higher exited states and the continuum should be dropped 
out by introducing the threshold sq as the lower bound of the integration. Since we are investigating 
the mixing, we insert all the three lowest states I/2) and I/3) into Eq.® and then we have: 

-ImWKs) = J2 mJi\fn){fn\mS{s-ml)+p'{s)e{s-s',), (16) 

n=l,2,3 

where i, j = q, s, g stand for the different currents (see Eq.©), n labels the state in the complete 
set, p^ {s) represents all the higher exited states and the continuum and Sq is the threshold for these 
higher states. 

Putting Eq. (fT6l) back into Eq.® and with the quark-hadron duality, we finally have the mo- 
ments as: 

= (0|Jq|/i)V-^ + (0|JJ/2)V-^ + (0|JJ/3)V-^ 

V^2e^™?" + ^'e"™^" + Vsle-'^'A {0\J^\Nf (17a) 



7^^^ = (o|Js|/i)V"^^ + (o|Js|/2)V™^- + (o|Js|/3)V"^^ 
= (0|Jg|/i)V"^^ + (0|Jg|/2)V™^ + (0|Jg|/3)V"^^ 

y^le^"?- + ^36-"*^^ + V.le-^^'A {0\J,\G)' (17c) 
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nf = (0|J,|A)(MJg|0)e-?^ + (0|J,|/2)(/2|Jg|0)e-i^ 
+ (0|J.|/3)(/3|Jg|0)e"™^ 
= (^VnV^ise-™?^ + ^21^236-"*^" + V^siV^sse-"^'^^ (0| JjiV) (G| JJO) (17d) 

= (0|Js|/i)(/i|Jg|0)e-^ + (0|Js|/2)(/2|Jg|0)e-™^ 
+ (0|Js|/3)(/3|^g|0)e-^ 
1^12^136-'"'^ + V22V23e-^"^ + Vs2Vs3e-"'"A (0| J, | ^) (Gl Jg|0) , (17e) 



where mi, m2 and rris are the masses of I/2) and I/3). In Eq.(fT7]). and the relationship of the 
physical and the un-physical states is involved in the calculations: such the concerned current only 
couples to the certain un-physical state with the right quantum number and flavor. For example, 
the current of the glueball cannot couple to the state of the light-quark, vice versa. For the physical 
state 1/1), 1/2) and I/3), we have: 

' {0\J,m = {0\J,\Va\N) = Va{0\J,\N) ; 

{0\Js\h) = {0\Js\V,2\S) = V,2{0\Js\S) ■ (18) 
[ (0| JJ/,) = {0\J,\V,s\G) = V,s{0\Ja\G) , 

where i = 1,2,3 for the three physical states. The un-physical states |A^), \S) and \G) directly 
couple to the certain currents, but do not correspond to any physical values. Thus we need to relate 
them to the physical states in terms via the moments in Eq. (fTTI) . Thus we are able to establish the 
equations for the ratios among the moments: 
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T^TA^.sf) _ rnVi3e-"^m^^^+^) + V-2iV-23e-^-m^^^^^^ + V-3iV-33e-"'^m^^^+^^ 

^0^) VuVi^e-^^mf + V^2iV23e-'^i"mi^ + V^iV^^e-^l^mf 

7^fc"+l(r, sf) _ V^^V^^e-^l-mf^^^ + V^22V^23e-'"^"m^^'+'^ + V^32^33e-'"3-m^ 
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(19) 
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FIG. 3: The Gray Area is the value of the Ihs. of the first tlir-ee equations in Eq.(fT9l):(a) soqq £ 
[3.1,3.7]GeV2; (b) sqss G [4.2, 4.8]GeV2;(c) sogg G [3.5, 4.1]GeV2. The Area between the Dashedline 
and the Dotted line is the value of the rhs. of the first thr^ee equations in Eq.([T9l): (a) [141,21 1 G [0.6, 0.8]; 
(b) |142,22| G [0.01,0.5]; (c) 1^13,231 G [0.5,0.8]. So the overlapping region is the proper parameter area for 
the mixing matrix V. 



Totally we have eight equations in Eq.(fT9l) and Eq.(fT4l) for determining the mixing matrix. Suppos- 
ing the matrix is real, there should be nine independent elements, but we only have eight equations, 
so that this equation group is not enough to directly determine the whole matrix. However, as we 
know, the matrix is unitary (as the matrix is real as assumed, it is an orthogonal matrix), thus we 
may gain an extra equation to fix all elements of the matrix. Namely, on the other hand, if we fix 
one element of the matrix V, in our work, for example, V23, then all other elements of the matrix 
V can be obtained by solving these eight equations. Sequently, let V23 run in the region [—1,1], 
the unitarity condition may help to eventually fix its value and the best fitting of the V is expected. 



IV. NUMERICAL RESULTS 

n Eq.(fT9l). it needs the values of the condensates and some other parameters as inputs. From 
[|33l1 . we set them as: 

rrig = O.OOSGeV, = 0.14GeV, 

mo = VaSGeV, {qq) = -0.24^GeV^ (20) 
(asG^) = 0.06GeV^ (^,05) = ml{qq)GeV^. 
The other parameters are related to the QCD sum rules: the Borel parameter r and the threshold 

of 

'^o'^' -^o^' -^o^' -^o^ ^iid ■^o'^ defined in Eq.([19l). By the general strategy, one should search for 
plateaus in the diagrams of the correlation versus the Borel parameter and the threshold sq. Only 
the parameters fall in a certain region, the plateaus can appear, namely within the plateaus the 
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FIG. 4: The Gray Area is the value of the Ihs. of the last two equations in Eq.([l9ll:(d) soqg G [3.6, 4.2]GeV^; 
(e) sosg G [4.2, 4.8]GeV^. The Area between the Dashed line and the Dotted line is the value of the rhs. of 
the last two equations in Eq.^: (d) Vn^2i G [-0.7,-0.6], V31 G [0.6,0.72], and V32 G [-0.72,-0.6]; 
(e) V21 G [0.1,0.3], V22 G [0.45,0.47], V31 G [0.6,0.7] and V32 G [-0.69, -0.60] 

results are not sensitive to the choice of Borel parameter and sq, then are trustworthy. In this 
work, there are six correlation functions in total, so we require all of them to have a common 
plateau region for the Borel parameter, where all the six moments are relatively independent of 
the Borel parameter. Obviously this condition is not easy to be satisfied. Once such a region is 
found, we would be able to conclude that the results based on the QCD sum rules make sense. The 
dependence of all six moments on the Borel parameter are presented in Fig -l3] and Fig-|4l And we 
can see obvious appearance of plateaus. 

We first have to check if in the parameter regions Eqs.(fT9^ have real solutions. We find that 
there are indeed. As we require the matrix V to be real, only a very narrow parameter space is 
available. The Fig-|3]and Fig-|4]show the values of the right-hand side (rhs) and the left-hand side 
(Ihs) of the equations in Eq.(fT9l). where the Fig-|3]is for the first three equations, and the Fig-|4]is 
for the last two equations. Taking the error tolerance into account, the lines would be widened into 
bands, in the Fig-|3]and Fig-|4l the region between the Dashed line and the Dotted line is for the rhs 
of the Eq.(fT9l) and the Gray one is for Ihs. It is clear that, only in the overlapping region, rhs and 
Ihs can be equal, and appearance of the overlapping region implies that a solution of the Eq.(fT9]) 
may exist. 

Searching for such an overlapping region in Fig-|3] and Fig-@] one needs to find a proper pa- 
rameter space. Eventually, we have found a satisfactory region where the best-fitted parameters 
are: the Borel parameter r G [1/1.8^, l/2.1^]GeV~^ and the five thresholds which must be close 
to s[J^ = 3.4GeV^ are s'^ = 4.5GeV^ sf = 3.8GeV^ = 3.9GeV^ and = 4.5GeVl At 
the same time, the allowed value ranges of the matrix elements Vij are also set. From Fig-|3l one 
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FIG. 5: The dependence of the matrix elements Vij on the Borel parameter r: in (a,b,c) the hne is for 

Vii^i2,i3, the dashed line is for V21, 22,23 and the Dotted line is for ^31^32,33 

notices that only as the matrix elements fall in the following regions: 

Vn E [-0.6, -0.8] V12 e [0.01, 0.5] V13 E [0.5, 0.8] 
V21 G [-0.6, -0.8] V22 e [0.3, 0.6] V^23 e [0.5, 0.8] 

(21) 

all the requirements are satisfied. It is also noted that due to the unitarity condition (fT4]| . V^i ( 
i = 1, 2, 3) depend on other elements Vu and V2i, thus their value-ranges would be uniquely 
determined (there might be a sign difference), once the others are fixed. 

Fig-|4]corresponds to the last two equations in Eq.(fT9]). and apparently overlapping regions exist 
when the matrix elements of V reside in the ranges (1211 . Moreover, for the last two equations in 
Eq.(fT9l). we set A; = 3. The reason is that, only when A; = or 3, the equations Eq.(fT9l) have real 
solutions. However, /c = is not proper since when k = 0, the Ihs does not appear in the plateau. 

We solve the equations Eq. (fT9l) together with the three equations in Eq. (fT4)) . Our strategy is to 
set V23 as a free parameter and let it run within a range. We find that only when V23 ~ —0.69, the 
matrix V is real and orthogonal. The numerical solution is given in Tab-U and the dependence of 
the matrix elements Vij on the Borel parameter r is shown in Fig-|5l 

Since such terms V^j exist in Eq.(fT9l). the solution may not be unique. As a matter of fact, 
we obtain eight independent groups of solutions. However, enforcing the unitary condition to the 
matrix V, we find that several groups are practically identical (i.e. they deviate from each other by 
just a common phase) and others must be dropped out because they do not satisfy the orthogonal 
condition. Finally only one group of solutions remains which is presented in the following table. 



1.85 1.90 1.95 2.00 2.05 2.10 
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1 /r (GeV) 
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TABLE I: The relationship of the elements of the V and r 

/l7^(GeV) 1.8 1.85 1.9 1.95 2.0 2.05 2.1 2.15 

Vii -0.71 -0.72 -0.73 -0.74 -0.74 -0.75 -0.75 -0.76 

V21 -0.72 -0.69 -0.65 -0.63 -0.60 -0.57 -0.55 -0.52 

^31 0.01 0.11 0.20 0.26 0.30 0.33 0.36 0.38 

V12 0.31 0.25 0.19 0.14 0.09 0.03 -0.06 -0.15 

V22 0.40 0.44 0.46 0.47 0.47 0.47 0.44 0.41 

V32 0.87 0.86 0.87 0.87 0.88 0.88 0.89 0.89 

^13 0.71 0.69 0.67 0.65 0.63 0.62 0.60 0.59 

V23 -0.69 -0.69 -0.69 -0.69 -0.69 -0.69 -0.69 -0.69 

V33 0.12 0.22 0.27 0.32 0.35 0.38 0.40 0.42 



Our numerical results show that r = l/2^GeV^ is the center of the common plateau, and the 
mixing matrix V is 



^ -0.74^°:°^ 0.09t°;| 0.631°;°^^ 

0.03 



V 



'-'■"U_o.08 -0.05 



-0.69* 



V 



U-'JU„o.19 



(22) 



0.881H^ 0.35«:?3V 



The numerical analysis indicates that the matrix elements Vu, V21, V22, V32 and V13 do not 
change much when the Borel parameter runs from l/1.8^GeV~^ to l/2.15^GeV~^, but it is also 
noted that the errors of V31, Vu and V33 are relatively larger. 

The ratio of the contribution of the perturbative part to the "Moments" TZ and the lowest state 
below the threshold is given in Fig-|6l 



V. CONCLUSION AND DISCUSSION 

From the FIG-16] we find that, within the range ofr G [1/2.152, l/1.82]GeV"^ the fraction 
of the perturbative part in the total contribution is over 60%. By the general principle of the 
QCD sum rules, after performing the Borel transformation, the perturbative contribution should 
dominate, and it is a criterion forjudging the reliability of the results. 60% is not too bad at all. 
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FIG. 6: The ratio of the contribution of the perturbative part to the "Moments" TZ and the lowest state under 
the threshold point 



For a comparison let us write down the mixing matrix given by Close et al. [|l8l] : 



c 



( 



\ 



-0.79 -0.13 0.60 
-0.62 0.37 -0.69 
0.14 0.91 0.39 



(23) 



In this work, we calculate the mixing of the | A^) , \S) and | G) to result in the physical resonances 
/o(1370), /o(1500) and /o(1710). The mixing matrix Eq. (|22l) which we obtained in the QCD sum 
rules is consistent with Vq lilSll which was achieved based on phenomenological studies. 

This work is based on the conjecture of Close and Kirk[|l8|l that only the mesons heavier than 
IGeV are mixtures of gg, ss and glueball G, because the lattice results indicate that the mass of 
O"*""*" is around 1.5 ~ 1.7GeV. This was also suggested by Narison et al. in their earlier papers [|35 - 
381]. 

With this picture we calculate the mixing off-diagonal correlators which result in the physical 
resonances /o(1370), /o(1500) and /o(1710). Narison and his collaborators computed the off- 



diagonal corre 



glueballil 



3 



ators in their pioneer work when they considered a mixing between meson and 



Today, thanks to the progress of experimental facilities and innovation of the data-analysis, 
many new resonances have been observed and data are updated. The available new data enable 
us to re-study the mixing effects, even though the basic techniques have been provided in those 
pioneer papers. That is the aim of this work. We are indeed very encouraged by the consistency 
between the numerical results obtained in terms of the QCD sum rules and that gained by the 
phenomenological research. It implies that the QCD sum rules are really a good approach for 
studying hadron physics even though certain uncertainties unavoidably exist. 
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Moreover, as we indicated above, the another two resonances /o(1790) and /o(1812) were 
observed and they also reside in the range of 1 to 2 GeV, therefore we do not have reason to ignore 
a possibility that all the five physical states /o(1370), /o(1500), /o(1710), /o(1790) and /o(1812) 
are mixtures of qq, ss, qqG and ssG and glueball of 0"*^^. But it would be much more difficult 
to calculate the mixing not only because then we have to deal with a five-dimensional matrix, 
but also the leading order of the perturbative part of the correlation function is two-loop feynman 
diagrams. But if it is the real physics, we need to carry out the calculations, and it will be the task 
of our next work. 
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